Chiral Kinetic Theory 
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We derive the non-equilibrium kinetic equation describing the motion of chiral massless particles 
in the regime where it can be considered classically. We show that the Berry monopole which appears 
at the origin of the momentum space due to level crossing is responsible for the chiral magnetic and 
vortical effects. 
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Introduction — The generation of non-dissipative cur- 
rents in a chiral (parity violating) system in response 
to an external magnetic field has attracted significant 
amount of interest recently. Such an effect, noted ear- 
lier in different contexts in Refs.[l|, has been recently 
proposed as an intriguing explanation for the charge- 
dependent correlations in heavy- ion collisions in Ref.[3|,|4( 
and termed chiral magnetic effect (CME). It has been 
also shown recently in Ref. [B| that hydrodynamics of chi- 
ral systems with anomaly requires the presence of such 
currents, as well as currents induced by the vorticity of 
flow - the chiral vortical effect (CVE), discovered earlier 
in a microscopic calculation in astrophysical context in 
Ref. @ and rediscovered recently in gauge-gravity duality 
calculations in Refs.@, H|- 

The interesting applications of these chiral transport 
effects involve highly non-equilibrium conditions, such as 
those arising in the early stages of the heavy-ion colli- 
sions, when the magnetic fields created by passing ions 
are still strong. However, derivations of these effects have 
been mostly done assuming thermal and chemical equi- 
librium. The aim of this paper is to address this short- 
coming of theory. 

A natural framework to study non-equilibrium condi- 
tions is a kinetic theory. As any useful theory, it has 
limitations, such as assumption of the classical motion 
between collisions and the weakness of the coupling. Nev- 
ertheless, a kinetic description would undoubtedly be an 
important step for our understanding of the chiral trans- 
port phenomena. 

Most of the ingredients of the approach presented here 
can be found in the literature on thephysics of geometric 
phases introduced by Berry in Ref. [9] . The relevant clas- 
sical equations of motion were introduced in Ref. [l(| (see 
also [11| for a review). The kinetic equation in the pres- 
ence of the Berry curvature have been studied, e.g., in 
Refs.fH 

[Hi ] . The most recent and closely relevant appli- 
cations include Refs. Unrelated to the above so 
far, but very important step towards kinetic description 
of the CME and CVE, was made recently in Ref. [161] . 

Putting these ingredients together we derive the de- 
sired non-equilibrium expressions for the CME and CVE 
which are, to the extent of our knowledge, new. The key 
observation of the present paper is that for Weyl fermions 
the Berry curvature, being the field of a monopole, leads 
directly to the chiral magnetic effect. We also point out 



that the chiral vortical effect can be similarly understood 
by simply replacing the Lorentz force due to the mag- 
netic field by the Coriolis force. In the following we shall 
present a reasonably self-contained derivation of these 
results using a formalism somewhat complementary to 
traditional approaches. This will allow us also to make 
connections to other field-theoretical concepts (such as 
abelian projection) more familiar in the particle theory 
context. 

Kinetic equation — Kinetic equation describes the mo- 
tion of particles in the regime where collisions are rare 
enough that motion between collisions is classical. In 
terms of the distribution function f(t,x,p) the equation 
reads 



df Of df df . 



(i) 



We think of a "cloud" of particles each of which follows its 
classical trajectory x(t), pit). As a result the distribution 
evolves with time in such a way that if one follows a local 
volume occupied by a set of particles along the trajectory, 
the number of particles in it can only be changed by 
collisions. 

Since collisions are not the focus of this paper we shall 
ignore collision integral C[f] in Eq. (QJ here. 

The CME is known to be closely related to the chiral 
anomaly @, @. On the other hand, it is clear from 
the description above that the number of particles in 
the phase space cannot change. How could a kinetic 
equation account for anomalous particle number non- 
conservation? In other words, how can classical equation 
account for quantum anomaly? As we shall see below, 
the answer is, in two words: Berry monopole. 

Path integral and U(2) gauge invariance — Consider 
the Hamiltonian for a Weyl particle: 



H = cr • p . 



(2) 



For each momentum p it represents a two-state system 
with energy gap 2\p\. 

It is more straightforward to obtain classical limit in 
the path integral formulation rather than in the canonical 
formulation of quantum mechanics usually employed to 
describe Berry connection. Consider the transition am- 
plitude between the two spin states i and /. By inserting 
the sums over complete sets of eigenstates of coordinates 
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and spin, \x,s), and momenta and spin, |p, A), the am- 
plitude can be written as a path integral 



</l 



T>xT>p V exp 



■p)dt 



(3) 



fi 



where we need to take a matrix element [. . .]fi of the 
path-ordered product of the matrices exp{— icr-pAt} over 
each path x(t), p(t) in the phase space. These matrices 
can be thought of as describing the rotation of the state 
of the particle in the spin space as it moves along. 

The massless particles we describe have only one he- 
licity state (the opposite helicity state corresponds to an 
antiparticle) . In order to consider classical motion of such 
a particle we need to diagonalize the matrix in the helic- 
ity basis and then apply the usual method of stationary 
phase to determine the classical trajectory. This diago- 
nalization can be done at each point on the trajectory 
using unitary matrix V p such that 



V^CT 



pV p = \p\a 3 . 



If the values of momenta at two neighboring points t\ and 
t2 are pi and pi , we insert identity matrices between the 
exponential factors in the following way: 

• • ■ V P2 V^ exp{-i<T ■ p 2 At}V P2 Vl 

x V Pl Vl exp{-z<T • Pl At}V Pl Vl . . . 
= ...V P2 e^{-i\p 2 \a 3 At}V^V Pl 

x e^p{-i\pi\a 3 At}Vl ... . (5) 

If the Ap = P2 — Pi is small, we can write the extra 
unitary rotation between the two neighboring exponents 
as 

V p2 V P i ~ exp(-ia p • Ap) , where a p = iV^V p V p (6) 

is a hermitian 2x2 matrix. 

Performing the above diagonalization along the whole 
trajectory and assembling the exponents into the path 
integral one obtains alternative expression for the ampli- 
tude in Eq. (0) 



(f\e 



xT^exp 



iH(t f -ti) 



V, 



Pf 



VxVp 



(p-x- \p\a 3 - a p • p)dt \ \£ 



(7) 



fi 



If we did not insist on diagonalizing the matrix er • p, we 
could have chosen an arbitrary U(2) rotation, say V p U p , 
instead of V p . This results in a local "gauge transforma- 
tion" of the "action" such that 

- \p\°3 -> -\p\U^a 3 U p , a p ->■ UlSipUp + iU p V p U p . 

(8) 



This gauge freedom corresponds to the free choice of the 
phase and spin quantization direction for the momentum 
states: |p, s) — > U p \p 7 s) along the trajectory. Clearly 
this choice only affects the expression for the amplitude 
in Eq. J7]), and not the value of the amplitude itself. We 
use this redundancy of description to choose the helicity 
basis at each p. This choice diagonalizes er-p and enables 
us to take the classical limit. 

Abelian projection and Berry monopole — Fixing this 
nonabelian U(2) gauge freedom by diagonalizing the 
Hamiltonian is mathematically similar to the abelian pro- 
jection introduced by 't Hooft in Ref.[17|. In the clas- 
sical regime the contribution of the transitions caused 
by the off-diagonal components of a p is negligible (in 
Ref. 17| the "non-abelian" part of the gauge field is non- 
propagating due to confinement). Suppressing these off- 
diagonal components, we still have a U(l)x U(l) gauge 
freedom corresponding to selecting arbitrarily the com- 
plex phases for the helicity eigenstates at each momen- 
tum. Focusing on helicity +1 we can denote the cor- 
responding diagonal component [a p ]n = a p . Then the 



(4) classical action for the helicity +1 particle becomes 



(p • x - \p\ - a p • p)dt . 



(9) 



The classical, or adiabatic, approximation will break 
down when the two eigenvalues of the Hamiltonian are 
degenerate, i.e., at p = 0. As we shall see, this point is 
the source of the effects of the quantum anomaly. 

As in the 't Hooft's original application of the abelian 
projection, even if the non-abelian field a p is a pure 
gauge, Eq. ([6]), the abelian component [a p ]n = a p is 
non-trivial. Finding the unitary matrix V p in Eq. (j4j and 
calculating a p in Eq. ©, one obtains the well-known re- 
sult that the corresponding abelian field a p is the field of 
a "monopole" TJl l]| at |p| = 0. Of course, the physical 
amplitude cannot depend on the gauge choice in Eq. © . 
We expect physical observables to depend only on the 
abelian field strength b — V p X a p . One finds 



b 



P 



2P 



whe 



re p : 



P_ 

\P\ 



(10) 



Equations of motion — Before we write the classical 
equations of motion, let us quantify the conditions of 
their applicability. The classical, or adiabatic, approxi- 
mation requires the off-diagonal components of a p • p in 
Eq. (J7J) to be small compared to the energy gap 2|p|. This 
means the forces, p, on the particle cannot be too strong. 
For example, if the particle moves in a magnetic field B: 
B < |p| 2 , where we used \a p \ - l/|p| (cf. Eq. (fTO)) ), 
This condition is obvious physically, since particles with 
momenta as low as the momenta on the lowest Landau 
orbit cannot be treated classically. 

It is easy to couple the classical particle described by 
the action in Eq. (0) to external electromagnetic field 
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given by scalar and vector potentials $ and A. By vari- 
ations of the resulting action 



Using notations E = \p\ and an overbar to denote av- 
erage over the unit sphere of directions of vector p we 
can write 



1=1 (p-x + A-x-$-\p\- a p -p)dt (11) 



j C ME = B f(p 



f(E,p)dE, (19) 



one obtains the desired equations of motion (cf. [HI ED]): This equation agrees with the result conjectured in 



x = p + p X b; 
p = E + x X B. 



(12) 
(13) 



Without the Berry flux b, these equations are familiar 
equations for the velocity of a massless particle and the 
Lorentz force. Without electromagnetic field the Berry 
curvature b drops out of the equations, because p = 0. 

Substituting Eq. (fTB")) into Eq. (TT2")) and solving for x 
one finds: 

VGx=p + E xb + B{p-b); (14) 
VGp = E + pX B + b(E-B). (15) 

Here G = (l+b-B) 2 is the determinant of the 6x6 matrix 
of coefficients in Eqs. (fT"2"|) . (|13|) for x and p. Substitut- 
ing into Eq. (fTJ) we can then obtain the desired kinetic 
equation for the distribution function of such particles in 
the phase space. 

Chiral magnetic effect — It is important to take ac- 
count of the fact that the invariant measure of the phase- 
space integration is given by \[G d 3 x d 3 pj (27r) 3 , see e.g., 
Ref . [20l| . In particular, one can check using equations 
of motion (Tl4l) . (fT5j) and Maxwell equations V • B = 0, 
Vx£ = dB I dt that this measure obeys Liouville equa- 
tion: 



2nE-BS 3 (p), (16) 



where the last term is due to the Berry monopole V p -6 = 
2ttS 3 (p), Eq. (TlUl) . The last term is the effect of the quan- 
tum anomaly which "injects" particle number violation 
into our otherwise classical description. It is notable that 
this term is localized at p = 0, where the classical descrip- 
tion must break down due to level crossing. 
The current density is given by 

i-frtf*, where j-jgfc. (17) 

Substituting Eq. ([13]) into Eq. (JT7J) we find 

j= [VGfx= f fp + Ex f fb + B f f(p.b). (18) 
J p J p J p J p 

The first term gives the usual current, while the second 
is the anomalous Hall current. Both vanish in a state 
with isotropic momentum distribution, such as equilib- 
rium state. The last term is the desired non-equilibrium 
expression of the CME. 



Ref. [16| for an isotropic distribution. In the case of the 



Fermi-Dirac distribution it reproduces the well-known re- 
sults (such as j'cme = f-B/i^ir) 2 at zero temperature). 

Chiral anomaly — To find the effect of the electro- 
magnetic anomaly we calculate the 4-divcrgence of the 
particle number current in Eq. (|17j) . It is illuminating 
to begin the discussion by introducing the 6+1 phase 
space current (p, px, pp), where p = \J~Gf obeys continu- 
ity equation with a source 

dp + d(px) + d(pp) 
dt 



dx 



dp 



2nE • B f 5 (p) 



(20) 



which follows from Eq. (TT]) (with C[f] =0 for simplicity) 
and Eq. (|16[) . Integrating over momentum p we obtain 



;E-Bf , 



(21) 



dn . 1 

L V • 7 = 

dt J 4tt 2 " 

where, as in Eq. (TT71) . (n,j) = J p {p, px) is the 3+1 space- 
time current density and fo is the value of the distribu- 
tion function / at p = 0. For the Fermi-Dirac distribu- 
tion at zero temperature and non-zero chemical potential 
/o = 1 and we reproduce the standard expression of the 
electromagnetic anomaly. 

Strictly speaking the above calculation is not com- 
pletely legitimate because we integrated over the whole 
momentum space, including the singular point p = 
where the classical description is not applicable. The 
way to think about this equation is to exclude the region 
\p\ < A around the singularity. The value of A must 
be large enough so that the classical description applies 
outside it (A > \/~B). 

Then, in the classical region \p\ > A, the 6+1 phase 
space current {p,px,pp) obeys continuity equation. I.e., 
the particles, in the absence of collisions, cannot be cre- 
ated or destroyed in the classical region. They can only 
enter or exit the region through the boundary of the re- 
gion at |p| = A. Integrating the continuity equation 
over the classical region \p\ > A and defining the 3+1 
current density in the classical region only (nA,jA) = 
^ p | >A (p, pi) we find that the non-conservation of the 
3+1 space-time current is matched by the momentum- 
space flux into the classical region through the boundary 
at |p| = A: 

dn A , „ , f dS A T ^ 



dt 1 V ' jA J (2tt)3 
where the flux density is given by 

J P = PP = (E + P X B) f + 2tt E • B f 



P 



47r|p| : 



(23) 
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The first term on the right-hand side is due to acceler- 
ation of the particles on the boundary |p| = A which 
moves them in or out of the classical region. This term 
gives a negligible contribution to the total flux in Eq. (j2"2")l 
if A is small enough that the variation of / over the 
boundary can be neglected. The total flux from the 
last term, however, tends to a finite limit when A — > 0, 
which is given by Eq. ([2"T]) . The origin of this net flux 
is the anomaly which operates, as is well-known, at the 
point of level crossing p = 0, lying inside the region 
\p\ < A, where the motion of particles must be treated 
fully quantum-mechanically. 

Chiral vortical effect — To describe chiral vortical ef- 
fect we need to realize that, unlike the external magnetic 
field B, which we can set directly, the vorticity w is a 
property of the flow of particles, which is indirectly con- 
trolled by external fields and initial conditions. More- 
over, the definition of vorticity involves hydrodynamic 
limit, which puts additional conditions on flow. How- 
ever, we can generalize the vorticity to non-equilibrium 
flows in the following way. We can decide to observe 
a given local fluid element in a co-moving frame, which 
will have to rotate with angular velocity ui with respect 
to the laboratory. The particles will experience addi- 
tional non-inertial forces in this frame. At this point we 
can generalize the problem to non-equilibrium by asking 
what additional currents such non-inertial forces induce. 

To linear order the only such force is the Coriolis force: 



p = 2\p\u X x + 0{ijj 2 



(24) 



(This classical result can be also verified by considering 
Weyl Hamiltonian in the rotating frame.) The effect of 
this force is the same as of a "magnetic field" B — > 2\p\u>. 
Making a corresponding substitution in Eq. (1151) we ar- 
rive at the following equation for the non-equilibrium 
generalization of the chiral vortical effect: 



i 



jcvE = w / 2\p\f(p-b) = 



Air 2 



f(E,p) 2EdE. 

(25) 

[16[ for isotropic 



This result is also in agreement with Ref. 
distribution, and reduces to j'cve = V?<^ I K llx T tor tlie 
well-known case of the Fermi-Dirac distribution at zero 
temperature. 

Conclusion — We presented kinetic description of the 
chiral magnetic and chiral vortical effects given by kinetic 
equation ([I]) with equations of motion (fl~2j) . fp~3|) . Al- 
though these equations are ubiquitous in the condensed 
matter literature on the effects of the Berry curvature, to 
our knowledge, their relationship to the chiral magnetic 
and chiral vortical effects has not been appreciated un- 
til now. The key observation that the Berry curvature 
for the Weyl Hamiltonian is sourced by a monopole at 
|p| = leads directly to the general non-equilibrium ex- 
pressions for the CME and CVE in Eqs. $F5§ and (|2"5)) 
which reproduce all known results in equilibrium. 



The presence of the monopole singularity in the mo- 
mentum space also provides a natural mechanism by 
which anomaly can change the particle number in an oth- 
erwise classical system. The classical description breaks 
down in the region surrounding the singularity at p = 
of the size of order of the typical momentum in the low- 
est Landau orbit. The net particle creation occurs by the 
purely quantum effect of anomaly (level crossing) inside 
this non-classical region. The net flux of the particles into 
the classical region is then given by Eq. (|22"j) . which can 
serve as a boundary condition for the kinetic equation in 
the classical region. 

It would be interesting to use the results obtained here 
to investigate the consequences of non-equilibrium for the 
chiral transport effects in heavy- ion collisions [H, 0, El . 
We leave this and other applications to further study. 
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